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Alexander Khaetskii^ , and Eugene Sukhorukov^ 

^ Department of Physics, University at Buffalo, SUNY, Buffalo, NY 14260-1500 

^ Department of Theoretical Physics, University of Geneva, 

24 quat Ernest Ansermet, CH-1211, Switzerland 

(Dated: February 24, 2013) 

We consider a 2D ballistic and quasi-ballistic structures with spin-orbit-related splitting of the 
electron spectrum. The ballistic region is attached to the leads with a voltage applied between 
them. We calculate the edge spin density which arises in the presence of a charge current through 
the structure. We solve the problem with the use of the method of scattering states and clarify the 
important role of the unitarity of scattering. In the case of a straight boundary it leads to exact 
cancellation of long-wavelength oscillations of the spin density. In general, however, the smooth spin 
oscillations with the spin precession length may arise, as it happens, e.g., for the wiggly boundary. 
Moreover, we show that the result crucially depends on the form of the spin-orbit Hamiltonian. 
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Currently, there is a great interest, both experimental 
and theoretical, to spin currents and spin accumulation in 
various mesoscopic semiconductor structures 1', '2'1. Both 
phenomena are due to spin-orbit (s-o) coupling and are 
of great importance for the future of spin electronics. 
The edge electron spin density accumulation, related to 
the Mott asymmetry in electron scattering off impurities, 
has been recently measured [3|. Moreover, the edge spin 
density in the two-dimensional (2D) hole system, which is 
due to the intrinsic mechanism i4!!| of the s-o interaction, 
has also been observed 5]. It is well known 2], that in 
the diffusive regime (and when a spin diffusion length is 
much larger than a mean free path), the spin density ap- 
pearing near the boundary is entirely determined by the 
spin flux coming from the bulk. For example, in the dif- 
fusive regime and in the case of the Rashba Hamiltonian, 
when the spin current in the bulk is zero, the spin density 
component perpendicular to the plane is zero everywhere 
down to the sample boundary [6]. 

In an opposite case, when the spin precession length 
is much shorter than the mean free path, the situation 
is much less investigated. An example of such a system 
is a mesoscopic structure with s-o-related splitting of the 
electron spectrum A^j, in the limit A^^Tp 3> 1, where Tp is 
the mean free time. Besides, in the presence of s-o inter- 
action, the boundary scattering itself is the source of gen- 
eration of the spin density. It is obvious, that the char- 
acteristic length near a boundary, where the spin density 
arises is the spin precession length, Lg = hvp/Aji, with 
vp being the Fermi velocity. This mechanism of the spin 
density generation is the subject of our paper. We show 
that various situations may arise, depending on the form 
of the s-o Hamiltonians. 

We start with a 2D system described by the Rashba 
Hamiltonian in the ballistic limit, where a mean free path 
is much larger than the sample sizes. The ballistic region 
is attached to the leads, and a voltage V applied between 
the leads causes a charge current through the structure. 



as shown in Fig.[T] Since the electric field is absent inside 
an ideal ballistic conductor, the edge spin polarization 
appears not as a result of the acceleration of electrons by 
an electric field, but rather due to the difference in pop- 
ulations of left-moving and right-moving electrons. The 
combined effect of boundary scattering and spin preces- 
sion leads to oscillations of the edge spin polarization. 
Note that there is no relation between the spin current 
in the bulk, which is zero in the considered case, and the 
edge spin accumulation. 

The problem of the spin density accumulation in a bal- 
listic system and for a straight boundary has been con- 
sidered earlier in Refs. [7|8| with the help of the Green's 
functions method. Surprisingly, in this case the long- 
wavelength oscillations of the spin density cancel, and 
the final result contains only Friedel-like oscillations with 
the momentum 2fci?. This effect may be interpreted as 
s-o splitting of the Friedel oscillations in the charge den- 
sity; two charge oscillations corresponding to spin-up and 
spin-down orientations get shifted with respect to each 
other in the presence of the s-o interaction. Therefore, 
strictly speaking, this phenomenon is different from a s-o- 
related accumulation of the spin density upon boundary 
scattering. Besides, the method used in Ref. [7| does not 
allow to understand the reason for the cancellation of 
long- wavelength oscillations of the spin density. 

We solve the problem of edge spin accumulation by 
using scattering theory, with scattering states coming 
from different leads of the structure and, therefore, hav- 
ing different occupations. The simplicity of the method 
allows us to gain an insight into the underlaying physics. 
We show that it is the unitarity of scattering that leads 
to the exact cancellation of long-wavelength oscillations 
of the spin density with the period L^ in the case of 
a straight boundary. It should be also mentioned that 
the observed behaviour is closely related to the effective 
one-dimensional character of scattering, arising from the 
translational invariance along the boundary. However, 
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FIG. 1: Left: Schematics of the boundary specular scattering 
in the presence of spin-orbit coupling. Plus and minus modes 
are shown for the same energy and the same wave vectors 
along the boundary. Right: Geometry of the system. Voltage 
V applied to the ideal leads causes a charge current through 
the ballistic region. 



the case of a straight boundary appears to be a rather 
exceptional one. In general, smooth spin oscillations with 
the spin precession length Lg arise, as it happens, for ex- 
ample, for the wiggly boundary, or for scattering off a 
circular impurity in a 2D electron system [9|, [lOJ . This is 
a consequence of the fact that in higher dimensions the 
conditions of the unitarity of scattering take a different 
form, as explained below. In all these situations, the spin 



density decays as a power law of the distance from the 
scatterer. 

Rashba s-o Hamiltonian in the bulk of a ballistic 2D 
electron system takes the following form 



^(p) -h^ ¥^' ^ p] 
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where n is the normal to the plane, a are the Pauli ma- 
trices, and p is the 2D momentum. The solutions of this 
Hamiltonian corresponding to the helicity values M = ± 
have the form exp(ipr/?i)xM(p), where r = x,y. The 
explicite form of the spinors and their eigenenergies is 
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with if being the angle between the momentum p and 
the positive direction of the a;-axis. 

We consider the semi-infinite system and choose the x- 
axis to be directed perpendicular to the boundary {x = 0) 
of the 2D system (see Fig.[T]). The wave functions, which 
obey zero boundary conditions at a; = 0, are obviously 
the scattering states, which constitute the complete set 
of the orthonormal functions. Two scattering states cor- 
responding to incident plus and minus modes with given 
wave vector along the boundary and the same energy are 
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Here, the wave vectors are defined as follows 



S acquire the following form: 



k^ = kl-kl kl=kl-kl, hk±=m{vFT^), (4) 
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(6) 



where p± = hk± are the momenta at the Fermi energy 
in the plus and minus modes. The angles ip, (pi may be 
expressed as sin((p) — ky/kj^ and sin(iy9i) ~ ky/k_ (see 
Fig. HI). 

From Eqs. ([2]) and ([3]), one finds the scattering ampli- 



tudes F^ and F^ 
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where Vxa = dci/dpx are the group velocities. For the 
Rashba model one has Vx,~/vx,+ — cos^pi/ cosip. 

The wave functions Eqs. ^ and ^ may now be used 
to calculate the average z-component of the spin as a 
function of coordinates: 
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One can check that the amplitudes FZ and F^ for the 
incident minus mode with the same ky and the same en- 
ergy are obtained from _F^ and Fl^ by replacing ip -^ ipi. 
Then, the components of the unitary scattering matrix 



where /i^(e, ky) is the Fermi distribution function, which 
takes either of two values, fF{€ — ^ — 6^/2) or fp{€ — 
ji -f eV/2), depending on the sign of ky. We find, that 
one may distinguish various contributions to {Sz{x)) with 
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FIG. 2: The original contour Fi along the real axis can be 
deformed into the part F2 going along the imaginary axis, 
and the part F3 going far from the origin. 



different oscillation periods, which originate from an in- 
terference of different terms in Eqs. (I2|) and ([3]). The 
smooth part of {Sz{x))s, which involves the interference 
of the outgoing waves [two last terms in Eqs. ^ and ^], 
reads: 



{Sz{x))s oc / dkydefpie^ky) 
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where 



^(X-(^i)l^.|x+M)e'^'-'^^" + ' 



A = 3+ ■ {s^y + Si ■ {szy . 



(8) 



Here we used the fact that the distribution function 
/^(e, ky), describing a particular lead, has the same value 
at given energy for plus and minus mode. Note, that the 
period of oscillations of the exponential factor e'C^"*^!)^ 
in Eq. (|S]) is of the order of the spin precession length. 
However, the term ([8]) vanishes, because the expression 
A is nothing but a non-diagonal component of the iden- 
tity matrix SS''. Thus, we obtain an interesting result, 
that the only reason for the cancellation of the long- wave 
length oscillations with the period Lg in {Sz{x)) is the 
unitarity of scattering. 

By taking into account in Eq. ^ the terms responsible 
for the interference between incoming and the outgoing 
waves [for example, between first and second terms in Eq. 
([2])], and adding the contribution from the evanescent 
modes [16|, we reproduce Eq. (16) of the Ref. [7|. It 
can be written in the form {Sz{x)) — {eV/8TT'^mVp)IinI, 
where 
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Note, that in the interval k^ < ky < k 



the quantity 

k has purely imaginary value, which corresponds to the 
evanescent modes. From this form of the presentation, we 
can immediately see that {Sz{x)) contains only 2kp com- 
ponent, while all the long- wavelength oscillations cancel 
exactly. Indeed, with the branch cut along the real ky 
axis between the points — fc__ and +k- (see Fig. [5]), the 



integrand function in / is an analytical function of the 
variable ky in the right half plane Keky > (for positive 
x). Since we need the imaginary part of /, the integra- 
tion is going alone the upper edge of the branch cut from 
up to +k- and then back along the lower edge of the 
branch cut. Because of the analyticity mentioned above, 
this integral is equal to the one taken along the imaginary 
axis of ky = in. Then, for x ^ Xp the latter integral is 
determined by small k <ti kp: 
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which gives for the spin density {Sz{x)) « 
{eV/2TT^VFx)cos{2mvFx)sm'^{amx/2), coinciding with 
the result of Ref. [7| . Therefore, the total spin per unit 
length along the boundary scales as /„ dx{Sz{x)) oc a^. 
Note, that the main contribution to this integral comes 
from small distances from the boundary, x c^ Xp- 

The cancellation of smooth spin density oscillations in 
case of the Rashba Hamiltonian and straight boundary 
occurs also in the quasi-ballistic situation: L 3> / 3> is, 
where L is the sample size, and / is the mean free path. 
In this case, the electric field in the bulk of the sample is 
finite. Therefore, the distribution functions for the plus 

and minus modes, /++(fc+) and / (k-), are determined 

by the electric field and by scattering off the impurities 
in the bulk of a system [ll|. The wave vectors fc+ and 
fc_ , shown in Fig. [l] correspond to a given energy and 
a given wave vector along the boundary. In the quasi- 
ballistic case considered here, these functions are equal, 
i.e. /++(fc+) = / (fc_), similar to a ballistic situa- 
tion. Under such a condition, the unitarity of scattering, 
see Eq.®, leads to the cancellation of smooth edge spin 
density oscillations, in contrast to what has been stated 



in the literature 13. Indeed, when the electric field is 



parallel to the boundary, the distribution functions in 

questions are /++(fc+) = f++(k+) s iny , and / (fc_) ~ 

J (fc_)sin(/3i (see Eq. (9) of Ref. [ll|). For the case of 

the Rashba Hamiltonian, the following relation has been 



|13j . Then, the ratio 
sin(^i = ky/ky — 1. 



obtained k+f--{k^) = fc_/++(fc+) 
is f++{k+)/f — (fc-) = fc+sin<^/fc_ 

Depending on the form of the s-o Hamiltonian, the 
unitarity may show up in totally different ways, leading, 
in general, to different patterns of the edge spin den- 
sity. In particular, boundary conditions play a crucial 
role. Let us consider 2D holes, described by the cubic (in 
2D momentum) s-o Hamiltonian, still keeping in mind 
the ballistic case and an abrupt straight boundary. The 
unitarity condition, i.e. the charge fiux conservation, ac- 
quires a different form. This may be seen already for 
the case of a normal incidence, where the Hamiltonian 
has the form p^/2m + a&xpl- In the case of the plus 
incident mode, applying previous zero boundary condi- 
tions (which sometimes are used for this Hamiltonian in a 
numerical treatment of the edge spin accumulation prob- 
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0, f: 



In other words, af- 




FIG. 3: Schematics of scattering of the plus incident mode 
by a wiggly boundary, x = VKsin(27rj//A). Apart from the 
main scattering channels (solid lines), there are additional 
scattering waves with the wave vectors along the boundary 
shifted by ±27r/A (dashed lines). 



ter elastic backscattering the helicity value changes sign, 
which is a consequence of the fact that a^ is conserved. 
Then, for the charge flux to be conserved, one needs the 
equality of the group velocities v+ and w_ , corresponding 
to the plus and minus modes at the same energy. How- 
ever, in contrast to the case of the Rashba Hamiltonian, 
where v+ — V- , those velocities are not equal for the 
cubic Hamiltonian: v^ — V- = a(p+ +p-)'^. The formal 
way to resolve the trouble is to note that the cubic Hamil- 
tonian has three solutions for a given energy, one of them 
corresponding to a momentum larger than kp (for small 
spin-orbit coupling). Thus, in general, the unitarity of 
scattering in slow plus and minus channels is violated, 
because some flux is carried away by a fast oscillating 
mode. Therefore, the corresponding coefficient A in Eq. 
dS]) does not vanish, and smooth spin density oscillations 
occur. However, for a physically relevant confining po- 
tential which grows to infinity in a smooth manner, the 
fast oscillating modes are not excited. Then, the uni- 
tarity of scattering into slow modes is restored, and the 
smooth spin density oscillations do not appear. 

Let us consider scattering off a wiggly boundary, shown 
in Fig. |3l for the case of the Rashba Hamiltonian. In this 
case the translational invariance is broken, therefore the 
condition of the unitarity of scattering takes a different 
form, as compared to the case of a straight boundary. 
As a result, the cancellation of the smooth spin density 
oscillations does not take place, leading to the total spin 



that is not small in the parameter a IJ, |15| . In order 
to demonstrate this effect, we consider the mathemati- 
cally simple case of the abrupt impenetrable boundary 
described by the equation: x = ({y) = Wsm{2ny/X). 
To the lowest order in W, the boundary condition reads: 



*(o,y) + C(y) 



dx 



(9) 



We are looking for the solution in the perturbative form 



,f,(o)^ 



,(1) 



^±{x,y) = \l/j_ (a;, 2/) -I- VP^ (cc,?/), where the zeroth or- 
der functions are given by Eqs. ^ and ([3]). The first 



order correction, proportional to W, is the superposi- 
tion of scattering waves with the wave vectors along the 
boundary shifted by ±27r/A (see Fig. [3]), and with the 
A:-vectors in the x-direction given by: 
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From now on, we assume Xp <^ L^ <t^ X. In ad- 
dition, in order to obtain an analytical expression for 
the spin density, we consider the case x/\/XXf ^ 1- 
In contrast to the case of the straight boundary, there 
are oscillations with three different periods: 2fci?- os- 
cillations, and the oscillations with two long periods. 



^ and Lg- Here ^ = 1/ \ k"^ — k^ is the new length 
scale. Under the conditions considered in the paper 
we obtain the set of inequalities Xp <C ^ ^ is, where 
Lg = l/(fc_ — /c+) — fi/ma is the spin precession length. 
If ky — )■ fc_|. (i.e., k — 7> 0), then ki tends to l/<^, which 
clarifies the physical meaning of ^. 

For the contribution of the long wavelength oscilla- 
tions, we obtain [Vt\: 
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,(10) 



where the last term is the contribution of the evanes- 
cent modes, and N\{x) is the Bessel function of the 
second kind. At the distances 2: ^ ^, we obtain the 
following dependence hong — —Ix^ jZS^ -I- (x/2i^)[7 + 
ln(a;/2Ls)]- In the opposite limit, a; 3> ^, we find 
hong = 2r7^i(^/^-) - ^ + ^^2£i£Zi). At even larger 
distances, a; S> is, one obtains smooth oscillations with 
the period of the order of Lg^ and the amplitude being 
proportional to ^fa. Note, that the total spin per unit 
length along the boundary is proportional to the integral 

/ dxhong{x) ~ {tt/2Ls) J^ dxNi{x/Ls) ~ J^ '^ dx/x ~ 

-ln(Ls/0 ~ hi(VAXF/0 - -(l/2)ln(A/AF), i.e., it is 
not small in s-o coupling, in contrast to the case of a 
straight boundary. 

In conclusion, we have considered the problem of the 
edge spin accumulation in mesoscopic structures with 
spin-orbit-related splitting of the energy spectrum, when 
the associated spin precession length is much smaller 
than the mean free path. In the presence of the charge 
current, the spin density develops oscillations near an 
edge in the direction transverse to the boundary. We 
have used the scattering states method to clarify the 
physics associated with this effect. The result crucially 
depends on the form of s-o Hamiltonian and the bound- 
ary conditions. The unitarity of scattering in the case 



of a straight boundary and Rashba Hamiltonian leads to 
the cancellation of long- wave length spin density oscilla- 
tions. On the contrary, the spin density in the case of 
wiggly boundary oscillates with a large period of the or- 
der of the spin precession length, quite similar to the case 
of electron scattering off antidots in 2D system with the 
Rashba Hamiltonian [10| . 

We acknowledge the financial support from the Pro- 
gram " Spintronics" of RAS, the Russian Foundation for 
Basic Research (Grant No. 07-02-00164-a), the Swiss Na- 
tional Science Foundation and SPINMET project (FP7- 
PEOPLE-2009-IRSES). 



[1] 



[2] 

[3] 
[4] 

[5] 
[6] 

[7] 

[8] 



H.A. Engel, E.I. Rashba, and B.I. Halperin, in Handbook 

of Magnetism and Advanced Magnetic Materials, ed. by 

H. Kronmiiller and S. Parkin, Vol.5 (John Wiley and 

Sons, New York, 2007). 

M.I. Dyakonov, and A.V. Khaetskii, in Spin Physics in 

Semiconductors, Springer Series in Solid- State Sciences, 

ed. by M.I. Dyakonov (Springer, Berlin, 2008). 

Y.K.Kato et al.. Science 306, 1910 (2004). 

S. Murakami et al.. Science 301, 1348 (2003); J. Sinova 

et al., PRL 92, 126603 (2004). 

J. Wunderlich et al., PRL 94, 047204 (2005). 

Ya. Tserkovnyak et al., PRB 76, 085319 (2007); O. 

Bleibaum, PRB 74, 113309 (2006). 

V.A. Zyuzin, P.G. Silvestrov, and E.G. Mishchenko, PRL 

99, 106601 (2007). 

G. Usaj and C.A. Balseiro, Europhys. Lett. 72, 621 

(2005); A. Reynoso et al., PRB 73, 115342 (2006). 



[9] Only in the effective ID situation the unitarity of scat- 
tering leads to the cancellation of long-wavelength os- 
cillations of the spin density. This is no longer the case 
for 2D systems. The problem for electrons, described by 
the Rashba Hamiltonian, scattered off a circular scatterer 
(an antidot) may be solved exactly in the presence of a 
charge current. The Sz-component of the spin density, 
oscillating with period Ls, appears around the antidot 

m- 

[10] A. Khaetskii, E. Sukhorukov, to be published. 

[11] A. Khaetsku, PRB 73, 115323 (2006). 

[12] E.B. Sonin, PRB 81, 113304 (2010). 

[13] Note the different meaning of quantities fc+ and k- used 
here, and quantiti es p + and p- used in [1J|. 

[14] Note, that Ref. [l^ considers the case of a straight 
boundary, which is smooth in the transverse direction, 
instead of abrupt one. If the width of the boundary is 
larger than Ls , two neighboring strips of the spin density 
of the opposite signs arise. We stress, that the physics 
discussed in Ref. [I^l is different from that we consider 
here. Since the boundary is still straight, the unitarity of 
scattering leads to the cancellation of smooth oscillations 
with the period La. As a result, the total spin, found in 
[151 ]. is small and proportional to a^. 

[15] P.G. Silvestrov, V.A. Zyuzin, and E.G. Mishchenko, PRL 
102, 196802 (2009). 

[16] When k+ < ky < fc_, the wave vector in x-direction is 
purely imaginary for the reflected plus mode. This mode 
is called evanescent. This situation takes place only for 
the incident minus mode, since only this mode comes 
from the leads under the indicated conditions. 

[17] For details, see A. Khaetsku, E. Sukhorukov, JETP Lett. 
92, 244 (2010). 



